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Superconductors host collective modes that can be manipulated with light. We 
show that strong terahertz light field can induce oscillations of the superconducting 
order parameter in NbN with twice the frequency of the terahertz field. The result 
can be captured as a collective precession of Anderson's pseudospins in ac driving 
fields. A resonance between the field and the Higgs amplitude mode of the 
superconductor then results in large terahertz third-harmonic generation. The 
method we present here paves a way toward nonlinear quantum optics in 
superconductors with driving the pseudospins collectively, and can be potentially 
extended to exotic superconductors for shedding light on the character of order 
parameters and their coupling to other degrees of freedom. 



Macroscopic quantum phenomena such as superconductivity and super- 
fluidity emerge in a variety of physical systems such as metals, liquid 
helium, ultracold atomic quantum gases, and neutron stars. One manifes- 
tation of the macroscopic quantum nature is the appearance of character- 
istic collective excitations. Indeed, phenomena associated with collective 
modes such as second sound and spin waves in condensates have been 
revealed in superfluid helium (1, 2) and in ultracold atomic gases (3, 4). 

Generally, collective modes in ordered phases arising from sponta- 
neous symmetry breaking are classified into (i) gapless phase modes 
(Nambu-Goldstone (NG) mode) and (ii) gapped amplitude modes (Higgs 
mode) (5-7). In charged-particle systems such as superconductors with 
long-range Coulomb interactions, the gapless NG mode becomes mas- 
sive, i.e., its energy is elevated to the plasma frequency as a result of the 
coupling to the gauge boson (photon field), which is referred to as the 
Anderson-Higgs mechanism (8, 9). The Higgs amplitude mode in super- 
conductors has been also studied theoretically (6, 70-75);because it is 
not accompanied by charge fluctuations, it does not couple directly to 
electromagnetic fields in the linear response regime. This is why the 
Higgs mode in conventional s-wave superconductors was observed only 
recently after a nonadiabatic excitation with a monocycle terahertz 
(THz) pulse (7(5); previous observations were in a special case where the 
superconductivity coexists with charge density wave that makes the 
Higgs mode Raman-active (77, 18). Hence many questions regarding the 
Higgs mode in superconductors remain unresolved: How does the mode 
couple to strong electromagnetic fields in nonlinear regime? Is it possi- 
ble to dynamically control the Higgs mode and therefore the supercon- 



ducting order parameter? 

Recent advances in the intense THz 
generation technique (19, 20) opens a 
new avenue for studying matter phases 
in non-equilibrium conditions. Ampli- 
tude- and phase-resolved spectroscopy 
using multi THz pulses has been real- 
ized (27), enabling the study of coher- 
ent transients in many-body systems in 
low energy range. Purpose of the pre- 
sent work is to explore coherent nonlin- 
ear interplay between collective mode 
in a superconductor and THz light field 
by investigating the real-time evolution 
of the order parameter under the driving 
field of a multi-cycle (as opposed to 
monocycle) THz pulse. 

In order to study evolutions on pi- 
cosecond time scale, we performed 
THz pump-THz probe spectroscopy 
(7(5, 22) (Fig. 1A). To generate an in- 
tense multi-cycle THz pulse as a coher- 
ent driving source, we first created an 
intense monocycle THz pulse by the 
tilted-pulse front method with a LiNb0 3 
crystal (19, 23). The monocycle pulse 
then goes through a band-pass filter to 
produce a narrow-spectrum multi-cycle 
pulse. Three band-pass filters are used 
to generate the different center frequen- 
cies at 0.3, 0.6, or 0.8 THz, respective- 
ly, with their power spectra displayed in 
Fig. IB. These photon energies are all 
below the superconducting gap of our 
NbN sample in the low temperature 
limit, which is 1.3 THz (Fig. 1C); this 
implies that the pump pulse does not 
generate quasiparticles (QPs) in one- 
photon processes at low temperatures. The sample is an s-wave super- 
conductor NbN thin film with 24-nm thickness grown on an MgO sub- 
strate (24) with r c =15 K. The ultrafast dynamics of the superconducting 
order parameter driven by the multi-cycle pump pulse is then probed 
through the transmittance of a monocycle THz pulse that enters the sam- 
ple collinearly with the pump pulse with a variable time delay. In gen- 
eral, we can detect the temporal waveform of the transmitted probe THz 
electric field E prohQ by varying the time delay of another optical gate 
pulse and using the electro-optic (EO) sampling method. In this experi- 
ment, we fix the timing of the optical gate pulse such that, in the absence 
of the pump, £ pro b e at this timing monotonically changes with tempera- 
ture, reflecting the change of the order parameter. Temporal evolution of 
the order parameter induced by the THz pump is sensitively monitored 
through the change of £ pro b e relative to its value in the absence of the 
pump as a function of the pump-probe delay time t pp (16, 22); we denote 
this change as ^ pro be- For details, see (25). Note that, in the present case, 
we could investigate the order parameter dynamics in the presence of 
coherently oscillating multi-cycle pump fields. The temporal waveform 
of the pump THz electric field E pump is displayed in Fig. ID for the cen- 
ter frequency of co=0.6 THz and the maximum electric field of 3.5 
kV/cm. Figure IE shows SE prohG with t pp at 7=14-15.5 K, at which co is 
greater than 2A(J). In this temperature range, the probe electric field 
gradually increases as a function of t pp to reach an asymptotic value, 
which indicates a reduction of the order parameter resulting from QP 
excitations (22). By contrast, at temperatures below 13 K where 2A(T) 
exceeds co (Fig. IF), the long-term reduction of the order parameter be- 
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comes less prominent as temperature decreases, because the QP excita- 
tion is suppressed. We can immediately notice that an oscillatory signal 
emerges with a frequency 1.2 THz (=2co) during the pump pulse irradia- 
tion, which indicates that the order parameter oscillates with twice the 
frequency of the driving field. 

We can physically capture the 2co oscillation of the order parameter 
in terms of the precession of Anderson's pseudospins (26). In the BCS 
ground state, two electrons with wavenumbers k and -k form a spin- 
singlet Cooper pair, with the BCS wavefunction given by 

i x f / Bcs)=nKi°°) k + ^i ii ) k )' (1) 

k 

where | oo) k ( 1 1 1) ) denotes unoccupied (occupied) k and -k. In Ander- 
son' s pseudospin formalism, the states h l\ and |00\ are represented by 

I / k I Ik 

up and down pseudospins, respectively (25), where the BCS ground state 
is thought of as a quantum superposition of up and down pseudospins 
with amplitudes v k and u k for each k, respectively (Fig. 2A). Note that 
the normal state has the pseudospins all up for \k\ < k ¥ and all down for 
\k\ > k ¥ at T=0 (Fig. 2B), where k ¥ is the Fermi wave number. Within this 
representation, the BCS Hamiltonian simply reads H BCS = 2^ b$ k » 

where a k = (a k x , o/, a k z ) is the Anderson's pseudospin (26) mapped onto 
the Bloch sphere (Fig. 2C), while b k =(-A\ -A", s k ) is a pseudomagnetic 
field acting on a h . Here e h is the band dispersion measured from the 
Fermi energy, A = A f +z"A"= U^^i^l + i^l) * s tne complex order pa- 
rameter, and U(>0) the pairing interaction. In equilibrium, each pseudo- 
spin is aligned along the pseudomagnetic field. The time evolution of the 
BCS state is then described by the Bloch equation for the pseudospins, 

— a lf = 2b k xa lf 

dt \(2) 
i.e., the time evolution of the BCS state is represented as the motion of 
the pseudospins in the pseudomagnetic field. For a spatially- 
homogeneous monochromatic electric field EQxp(icot) irradiated onto a 
superconductor, the z component of b k , in the nonlinear response regime, 
becomes (25) 



-Z 



dtdk 



-4(04(0 + 0(A 4 ) 



-tI 5s e,Ej exp( ' 2 ^ + ° {A ^' 

2 Tfdtdk, a> 



(3) 



where A(t) is the vector potential representing the electric field. We can 
see that the leading term in the energy variation is ~A(if, which is intui- 
tively because electrons with charge -e hybridize with holes with charge 
+e in the condensed pair, leading to a nonlinear coupling between light 
and the condensate, and to the pseudospin precession with the frequency 
2co. The coherent collective precession of the pseudospins (Fig. 2D) 
manifests itself macroscopically as the order parameter oscillation, and 
the change of the order parameter in turn affects the pseudomagnetic 
field. We calculated the time evolution of the order parameter self- 
consistently by numerically solving the Bloch equation for the multi- 
cycle pulse (see (25) for details). The simulation indeed exhibits the 
order parameter oscillation with twice the frequency of the external elec- 
tric field (Fig. 2E). 

An analytic solution for the linearized Bloch equation can in fact be 
obtained (27), where the temporal variation of the order parameter am- 
plitude [A(t)= A+8A(t)] turns out to behave around 2co=2A as 

SA(t) oc — cos(2## - <p) 

,(4) 



\2co-2A\ 



with ([) being a phase shift that depends on co. The divergence of the am- 
plitude at 2co=2A can be interpreted as a resonance between the induced 
pseudospin precession with frequency 2co and the collective amplitude 



mode of the order parameter, namely the Higgs mode (16), with fre- 
quency 2A. We can then relate the nonlinear current density y" NL in- 
duced by the external ac field with the change in the order parameter 
8A(0 via (25) as 

h L (t)K^A(t)SA(t) 

U ,(5) 

which takes a form of the London equation. This enables us to regard the 
nonlinear current, which reflects the dynamics of the order parameter, as 
a part of the supercurrent. 

Because the coherent interaction between the superconductor and the 
THz electromagnetic radiation results in the nonlinear (2co, 4co, 6co...) 
oscillations of the order parameter (Eq. 3), the nonlinear current in Eq. 5 
should accommodate higher odd-order harmonics in the transmitted 
pump THz pulse. To confirm this, we have performed a nonlinear trans- 
mission experiments for the pump THz pulse (without the probe pulse) 
(Fig. 3A). Figure 3B shows the waveforms of the transmitted pump THz 
pulse above (15.5 K) and below (10 K) r c =15 K; the waveform of the 
transmitted pulse below T c is considerably distorted. The power spectra 
of the transmitted pump THz pulse are shown in Fig. 3C on a logarith- 
mic scale and in Fig. 3D on an expanded linear scale at various tempera- 
tures. Below 7M3 K, a prominent peak appears around 1.8 THz, which 
indeed coincides with 3 go. The intensity at 3 co as a function of the pump 
electric field strength, depicted in Fig, 3E on a log-log scale, obeys 
l^pumpl 6 dependence, endorsing that the signal arises from the third har- 
monic generation (THG). The THG intensity at 10 K normalized by that 
of the incident pump pulse reaches 8x1 0~ 5 , which is high for a film with 
only 24-nm thickness and 3.5-kV/cm peak electric field (28). We can 
increase the interaction length up to about 0.2 um, the penetration depth 
of the sample at 0.6 THz (29), which would result in even higher conver- 
sion efficiency. A shift of THG peak energy with temperature is dis- 
cerned in Fig. 3D, which is attributed to the softening of the Higgs mode 
(2A(T)) toward T c . The THG signal disappears before the softening 
completes because the resonant enhancement is rapidly suppressed when 
A(T) moves out of the narrow bandwidth of the incident pump field. 

The 2 co oscillation of the order parameter and the THG are also ob- 
served for co=0.3 and 0.8 THz pumping. Figure 4B summarizes the tem- 
perature dependence of the THG intensities for all of co=0.3, 0.6, and 0.8 
THz. For co=0.3 and 0.6 THz, the THG signal peaks at 13.5 K and 10 K, 
respectively, whereas the THG signal for co=0.8 THz monotonically 
increases with decreasing temperature. Comparing the temperature de- 
pendence of the order parameter 2A(T) (Fig. 4A) with twice the pump 
frequency 2co (=0.6, 1.2, and 1.6 THz), one can deduce that the peak in 
the THG does fall upon 2co=2A(7). The THG intensity in Eq. 5 depends 
on the change of the order parameter amplitude, which is resonantly 
enhanced when 2co approaches the inherent Higgs amplitude mode 
2A(T). Indeed, the temperature dependence of the THG intensity calcu- 
lated with Eq. 5 shown in Fig. 4C agrees qualitatively with experiment 
in Fig. 4B. We conclude that the resonance of the Anderson's pseudo- 
spin precession in the superconductor is achieved by irradiation of THz 
pump, which results in large THG. We can note that the theoretical re- 
sults in Fig. 4C exhibit sharp resonance peaks, which results from the 
lifetime of the Higgs mode assumed to be infinite (i.e., power-law de- 
cay) within the BCS approximation (10, 12). In contrast, the observed 
resonance widths in Fig. 4B are finite, which may be caused by decaying 
channels for the Higgs mode and the finite spectral width of the pump 
pulse (Fig. IB). There are in fact various possible decay processes, in- 
cluding scattering with single-particle excitations, impurities, phonons, 
or low-frequency NG mode that emerges near T c (30), for which system- 
atic studies are desirable (31). 

We finally note that superconductors are known to exhibit highly 
nonlinear responses near the critical field or temperature, giving rise to 
nonlinear I-V characteristics and higher-order harmonics in transport 



Sciencexpress / http://www.sciencemag.org/content/early/recent / 10 July 2014 /Page 2/ 10.1 126/science. 1254697 



measurements with a frequency range from a few hertz to microwave 
(32-34). By contrast, the large nonlinear optical effect revealed here 
originates from resonance of ac fields to the collective amplitude mode 
of the order parameter, which leads to the strong THG emission in THz 
frequency range. 

The time-resolved observation of the THz higher-order harmonics 
will provide a unique avenue for probing ultrafast dynamics of the order 
parameter in out-of-equilibrium superconductors. It is highly intriguing 
to explore the quantum trajectories of the pseudospins on Bloch sphere 
in the non-perturbative light-matter interaction regime with much higher 
THz fields, which would result in a dynamics of superconducting order 
parameter not attained in conventional regimes. The present scheme 
utilizing the nonlinear coupling between pseudospins and light can be 
also extended to unconventional superconductors such as the cuprate or 
iron-pnictide, which would provide new insight about the high- r c super- 
conductivity and the interplay between the superconducting phase and 
other coexisting/competing orders. 
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Fig. 1. THz pump-THz probe spectroscopy. (A) Schematic experimental setup for the THz pump-THz probe spectroscopy, 
where BPF is a metal-mesh band-pass filter and WGP a wire-grid polarizer. (B) Power spectra of the pump THz pulse with 
the center frequencies of co=0.3, 0.6, and 0.8 THz. (C) Temperature dependence of the superconducting gap energy 2A of 
the NbN sample evaluated from the optical conductivity spectra based on the Mattis-Bardeen model (35). Horizontal lines 
indicate the center frequencies of the pump pulse. (D) Waveform of the pump THz electric field E pump with the center 
frequency of co=0.6 THz, with the squared |E pum p| 2 also shown. (E) The change in the probe THz electric field 5E pr0 be as a 
function of the pump-probe delay time £ pp in the temperature range 2A(7)<co. Increase of <5E pr0 be corresponds to a reduction of 
the order parameter. (F) 5E pr0 be against £ pp in the temperature range 2A(7)>co. 



Sciencexpress / http://www.sciencemag.org/content/early/recent / 10 July 2014 /Page 4/ 10.1 126/science. 1254697 




Time t 



Fig. 2. Anderson's pseudospin model and simulation with Bloch equation. (A, B) Schematics of the electron distribution 
represented by Anderson's pseudospins for the normal state at 7=0 and for the BCS state, respectively. (C) The pseudospins 
mapped on the Bloch sphere. (D) A schematic picture of the pseudospin precession. (E) Simulation of the Bloch equation 
showing the temporal evolution of the order parameter in an electric field (for details see (25)). 
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Fig. 3. THG in transmission spectroscopy. (A) A schematic of the nonlinear THz transmission experiment. (B, C) 
Waveforms and power spectra of the transmitted pump THz pulses below (10 K) and above (15.5 K) 7" c =1 5 K, respectively. 
(D) Power spectra of the transmitted pump THz pulse at various temperatures. (E) THG intensity as a function of the pump 
THz field strength. 
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Fig. 4. Temperature dependence of the THG intensity. (A) Temperature dependence of the order parameter 2A(7~) 
compared with twice the pump frequencies 2co (horizontal lines). (B) Measured temperature dependence of the THG 
intensities at co=0.3, 0.6 and 0.8 THz. (C) Calculated THG intensities as a function of temperature obtained by numerically 
solving the Bloch equation. 
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